*Uttermark Paper II Dissertation Project*
*UTTERMARK, MATTHEW J.*
*Edited 10/09/2017*



********INSTALL THE FOLLOWING IF NEEDED********
* tufte
* fsum
* estout
* tabout
* mplotoffset
* revrs
* cem


/* SET DIRECTORY.*/
**************
cd "C:/Users/muttermark laptop/Documents/stata/intersection/PAR_Replication"
**************
capture log close //  close any previous logs you had been working on
log using par_replication_062723, replace text //  log all work
version 16.1 // set version
clear all // clear memory

/* Load data set*/
use par_asap.dta

set scheme tufte


***Gen Gender Race Variable ***

gen intersection=0 if k_4a==1 & k_3a==1
replace intersection=1 if k_4a==1 & k_3a==2
replace intersection=2 if  k_4a==0 & k_3a==1
replace intersection=3 if  k_4a==0 & k_3a==2

label define mintersection 0"White Man"  1"White Woman" 2"Man of Color" 3"Woman of Color"

label value intersection mintersection


rename a_2a civil_service 
label var civil_service "Civil Servant"
rename a_8a weekly_hours
label var weekly_hours "Average Weekley Hours Worked"
rename j_1a years_employ_state
label var years_employ_state  "Years Employed in State Gov."
rename j_1b years_employ_agency
label var years_employ_agency  "Years Employed in Agency"
rename j_1c years_employ_position 
label var years_employ_position  "Years Employed in Position"
gen age_2 = k_2a*k_2a
label var age_2 "Age-squared"
gen age = sqrt(age_2)
label var age "Age"
gen pid5= 1 if k_8b==2
replace pid5= 2 if k_8b==4
replace pid5= 3 if k_8b==5
replace pid5= 4 if k_8b==3
replace pid5= 5 if k_8b==1
label var pid5 "Party ID"
label define mpid5 1"Republican"  2"Lean Republican" 3"Independant" 4"Lean Democratic" 5"Democratic"

label value pid5 mpid5


rename k_16a edu
label var edu "Education"

rename a_3b agency_size
label var agency_size "Total Agency Employees"

rename a_4b agency_budget 
gen log_agency_budget = ln(1+agency_budget) 
label var log_agency_budget "ln(Agency Budget, $2018)"

revrs e_1a 
revrs e_1b 
revrs e_1c 
revrs e_1d 

label var reve_1a "Gov."
label var reve_1b "Gov. Staff"
label var reve_1c "Legis."
label var reve_1d "Legis. Staff"

fsum reve_1a reve_1b reve_1c reve_1d  intersection civil_service weekly_hours age age_2 edu years_employ_state years_employ_agency years_employ_position pid5 agency_size log_agency_budget, stats(n mean sd min median max) catvar(  intersection)

*for table alt 1
*drop if state== 27 | state==34 | state==45 
*for table alt 2
*drop if state==11

**** CODE USED TO PRODUCE TABLE 1*******

fsum reve_1a reve_1b reve_1c reve_1d  d_15a d_16a d_20a d_21a intersection civil_service weekly_hours age age_2 edu years_employ_state years_employ_agency years_employ_position pid5 agency_size log_agency_budget, uselabel stats(min median max) catvar(reve_1a reve_1b reve_1c reve_1d d_15a d_16a d_20a d_21a intersection civil_service edu) 

**** TABLE EXPORETED AND FORMATTED IN EXCEL ****

***** CODE USED TO PRODUCE TABLE 2 ******

ologit reve_1a i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto m1

ologit reve_1b i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto m2

ologit reve_1c i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto m3

ologit reve_1d i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto m4

 esttab m1 m2 m3 m4 using Table_2.rtf ,starlevels(+ 0.10 * 0.05 ** 0.01) cells(b(star fmt(3)) se(par fmt(3))) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Gov." "Gov. Staff" "Legis." "Legis. Staff" )


***** CODE USED TO PRODUCE TABLE 3 ******
est restore m1
margins, dydx(3.intersection)

est restore m2
margins, dydx(3.intersection)

est restore m3
margins, dydx(3.intersection)

est restore m4
margins, dydx(3.intersection)

**** RESULTS EXPORTED AND FORMATTED IN EXCEL ****

**** CODE USED TO PRODUCE TABLE 3 ****

ologit d_15a i.intersection i.reve_1a i.reve_1b i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year , r 
est sto m5

ologit d_16a i.intersection i.reve_1c  i.reve_1d i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year , r 
est sto m6

 
ologit d_20a i.intersection i.reve_1a i.reve_1b  i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year , r  
est sto m7

ologit d_21a i.intersection i.reve_1c i.reve_1d i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year , r 
est sto m8 
 
 
  esttab m5 m6 m7 m8 using Table_3.rtf ,starlevels(+ 0.10 * 0.05 ** 0.01) cells(b(star fmt(3)) se(par fmt(3))) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Gov. Policy" "Legis. Policy" "Gov. Regs" "Legis. Regs")

***** CODE USED TO PRODUCE FIGURE 1 PANEL A ******

est restore m5
margins, dydx(intersection)
mplotoffset, recast(scatter) title(Governor) xtitle(Influence on Major Policy Decisions) xlabel(2 "White Woman" 3 "Man of Color" 4 "Woman of Color", labsize(small)) legend(rows(1) stack size(small)) legend(order ( 5 "None" 6 "Slight" 7 "Moderate" 8 "High")) aspect(0.9)

graph save figure_1a
 

***** CODE USED TO PRODUCE FIGURE 1 PANEL B ******

est restore m6
margins, dydx(intersection)
mplotoffset, recast(scatter) title(Legislature) xtitle(Influence on Major Policy Decisions) xlabel(2 "White Woman" 3 "Man of Color" 4 "Woman of Color", labsize(small)) legend(rows(1) stack size(small)) legend(order ( 5 "None" 6 "Slight" 7 "Moderate" 8 "High")) aspect(0.9)

graph save figure_1b

***** CODE USED TO PRODUCE FIGURE 2 PANEL A ******

est restore m7
margins, dydx(intersection)
mplotoffset, recast(scatter) title(Governor) xtitle(Influence on Rules and Regulations) xlabel(2 "White Woman" 3 "Man of Color" 4 "Woman of Color", labsize(small)) legend(rows(1) stack size(small)) legend(order ( 5 "None" 6 "Slight" 7 "Moderate" 8 "High")) aspect(0.9)

graph save figure_2a


***** CODE USED TO PRODUCE FIGURE 2 PANEL B ******

est restore m8
margins, dydx(intersection)
mplotoffset, recast(scatter) title(Legislature) xtitle(Influence on Rules and Regulations) xlabel(2 "White Woman" 3 "Man of Color" 4 "Woman of Color", labsize(small)) legend(rows(1) stack size(small)) legend(order ( 5 "None" 6 "Slight" 7 "Moderate" 8 "High")) aspect(0.9)

graph save figure_2b

***** TO CALCULATE TABLE NOTES FOR EACH GRAPH RUN THE FOLLOWING ****
*** 1A ***
est restore m5
test i0.intersection=i3.intersection
test i1.intersection=i3.intersection
test i2.intersection=i3.intersection
*** 1B ***
est restore m6
test i0.intersection=i3.intersection
test i1.intersection=i3.intersection
test i2.intersection=i3.intersection
*** 2A ***
est restore m7
test i0.intersection=i3.intersection
test i1.intersection=i3.intersection
test i2.intersection=i3.intersection
*** 2B ***
est restore m8
test i0.intersection=i3.intersection
test i1.intersection=i3.intersection
test i2.intersection=i3.intersection

*** RESULTS EXPORTED AND FORMATTED INTO FIGURES 1 & 2 ***

************************************************
*** REPLICATION CODE FOR SUPPLIMENTAL TABLES ***
************************************************

***** CODE USED TO PRODUCE APPENDIX TABLE B1 ******

ologit reve_1a b1.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto a1

ologit reve_1b b1.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto a2

ologit reve_1c b1.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto a3

ologit reve_1d b1.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto a4


est restore a1
margins, dydx(3.intersection)
est restore a2
margins, dydx(3.intersection)
est restore a3
margins, dydx(3.intersection)
est restore a4
margins, dydx(3.intersection)

**** RESULTS EXPORTED AND FORMATTED IN EXCEL ****

***** CODE USED TO PRODUCE APPENDIX TABLE B2 ******
ologit reve_1a b2.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto a5

ologit reve_1b b2.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto a6

ologit reve_1c b2.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto a7

ologit reve_1d b2.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto a8


est restore a5
margins, dydx(3.intersection)
est restore a6
margins, dydx(3.intersection)
est restore a7
margins, dydx(3.intersection)
est restore a8
margins, dydx(3.intersection)

**** RESULTS EXPORTED AND FORMATTED IN EXCEL ****

***** CODE USED TO PRODUCE APPENDIX TABLE B3 ******


drop if state==27 | state==34 | state==45

ologit reve_1a i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year  , r 
est sto b1
ologit reve_1b i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto b2
ologit reve_1c i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto b3
ologit reve_1d i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , r 
est sto b4
ologit d_15a i.intersection i.reve_1a i.reve_1b i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year , r 
est sto b5
ologit d_16a i.intersection i.reve_1c  i.reve_1d i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year , r 
est sto b6
ologit d_20a i.intersection i.reve_1a i.reve_1b  i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year , r  
est sto b7
ologit d_21a i.intersection i.reve_1c i.reve_1d i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year , r 
est sto b8 
 
 esttab b1 b2 b3 b4 b5 b6 b7 b8 using Table_B3.rtf ,starlevels(+ 0.10 * 0.05 ** 0.01) cells(b(star fmt(3)) se(par fmt(3))) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Gov. Policy" "Legis. Policy" "Gov. Regs" "Legis. Regs" "Gov. Policy" "Legis. Policy" "Gov. Regs" "Legis. Regs") 

clear all
***** RE-RUN PREAMBLE ***** 
use par_asap.dta

gen intersection=0 if k_4a==1 & k_3a==1
replace intersection=1 if k_4a==1 & k_3a==2
replace intersection=2 if  k_4a==0 & k_3a==1
replace intersection=3 if  k_4a==0 & k_3a==2

label define mintersection 0"White Man"  1"White Woman" 2"Man of Color" 3"Woman of Color"

label value intersection mintersection


rename a_2a civil_service 
label var civil_service "Civil Servant"
rename a_8a weekly_hours
label var weekly_hours "Average Weekley Hours Worked"
rename j_1a years_employ_state
label var years_employ_state  "Years Employed in State Gov."
rename j_1b years_employ_agency
label var years_employ_agency  "Years Employed in Agency"
rename j_1c years_employ_position 
label var years_employ_position  "Years Employed in Position"
gen age_2 = k_2a*k_2a
label var age_2 "Age-squared"
gen age = sqrt(age_2)
label var age "Age"
gen pid5= 1 if k_8b==2
replace pid5= 2 if k_8b==4
replace pid5= 3 if k_8b==5
replace pid5= 4 if k_8b==3
replace pid5= 5 if k_8b==1
label var pid5 "Party ID"
label define mpid5 1"Republican"  2"Lean Republican" 3"Independant" 4"Lean Democratic" 5"Democratic"

label value pid5 mpid5


rename k_16a edu
label var edu "Education"

rename a_3b agency_size
label var agency_size "Total Agency Employees"

rename a_4b agency_budget 
gen log_agency_budget = ln(1+agency_budget) 
label var log_agency_budget "ln(Agency Budget, $2018)"

revrs e_1a 
revrs e_1b 
revrs e_1c 
revrs e_1d 

label var reve_1a "Gov."
label var reve_1b "Gov. Staff"
label var reve_1c "Legis."
label var reve_1d "Legis. Staff"

***** CODE USED TO PRODUCE APPENDIX TABLE B4 ******
ologit reve_1a i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year if state!= 11 , r 
est sto c1
ologit reve_1b i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year if state!= 11 , r 
est sto c2
ologit reve_1c i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year if state!= 11, r 
est sto c3
ologit reve_1d i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year if state!= 11 , r 
est sto c4
ologit d_15a i.intersection i.reve_1a i.reve_1b i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year if state!= 11 , r 
est sto c5
ologit d_16a i.intersection i.reve_1c  i.reve_1d i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year if state!= 11 , r 
est sto c6
ologit d_20a i.intersection i.reve_1a i.reve_1b  i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year if state!= 11 , r  
est sto c7
ologit d_21a i.intersection i.reve_1c i.reve_1d i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year if state!= 11 , r 
est sto c8 

 esttab c1 c2 c3 c4 c5 c6 c7 c8 using Table_B4.rtf ,starlevels(+ 0.10 * 0.05 ** 0.01) cells(b(star fmt(3)) se(par fmt(3))) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Gov. Policy" "Legis. Policy" "Gov. Regs" "Legis. Regs" "Gov. Policy" "Legis. Policy" "Gov. Regs" "Legis. Regs") 


***** CODE USED TO PRODUCE APPENDIX TABLE B5 ******
gen treated =1 if intersection==3
replace treated=0 if intersection==0 

drop if treated==. 
cem   civil_service  age(#6) age_2(#6) edu(#0) funcat(#0), treatment(treated)


ologit reve_1a i.treated weekly_hours years_employ_state b3.pid5 agency_size log_agency_budget years_employ_agency years_employ_position inst6017_nom i.state i.year [iweight=cem_weights] , r
est sto d1
ologit reve_1b i.treated weekly_hours years_employ_state b3.pid5 agency_size log_agency_budget years_employ_agency years_employ_position inst6017_nom i.state i.year [iweight=cem_weights] , r
est sto d2
ologit reve_1c i.treated weekly_hours years_employ_state b3.pid5 agency_size log_agency_budget years_employ_agency years_employ_position inst6017_nom i.state i.year [iweight=cem_weights] , r
est sto d3
ologit reve_1d i.treated weekly_hours years_employ_state b3.pid5 agency_size log_agency_budget years_employ_agency years_employ_position inst6017_nom i.state i.year [iweight=cem_weights] , r
est sto d4
ologit d_15a i.treated weekly_hours years_employ_state b3.pid5 agency_size log_agency_budget years_employ_agency years_employ_position inst6017_nom i.state i.year [iweight=cem_weights] , r
est sto d5
ologit d_16a i.treated weekly_hours years_employ_state b3.pid5 agency_size log_agency_budget years_employ_agency years_employ_position inst6017_nom i.state i.year [iweight=cem_weights] , r
est sto d6
ologit d_20a i.treated weekly_hours years_employ_state b3.pid5 agency_size log_agency_budget years_employ_agency years_employ_position inst6017_nom i.state i.year [iweight=cem_weights] , r
est sto d7
ologit d_21a i.treated weekly_hours years_employ_state b3.pid5 agency_size log_agency_budget years_employ_agency years_employ_position inst6017_nom i.state i.year [iweight=cem_weights] , r
est sto d8

 esttab d1 d2 d3 d4 d5 d6 d7 d8 using Table_B5.rtf ,starlevels(+ 0.10 * 0.05 ** 0.01) cells(b(star fmt(3)) se(par fmt(3))) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Gov. Policy" "Legis. Policy" "Gov. Regs" "Legis. Regs" "Gov. Policy" "Legis. Policy" "Gov. Regs" "Legis. Regs") 


clear all
***** RE-RUN PREAMBLE ***** 
use par_asap.dta

gen intersection=0 if k_4a==1 & k_3a==1
replace intersection=1 if k_4a==1 & k_3a==2
replace intersection=2 if  k_4a==0 & k_3a==1
replace intersection=3 if  k_4a==0 & k_3a==2

label define mintersection 0"White Man"  1"White Woman" 2"Man of Color" 3"Woman of Color"

label value intersection mintersection


rename a_2a civil_service 
label var civil_service "Civil Servant"
rename a_8a weekly_hours
label var weekly_hours "Average Weekley Hours Worked"
rename j_1a years_employ_state
label var years_employ_state  "Years Employed in State Gov."
rename j_1b years_employ_agency
label var years_employ_agency  "Years Employed in Agency"
rename j_1c years_employ_position 
label var years_employ_position  "Years Employed in Position"
gen age_2 = k_2a*k_2a
label var age_2 "Age-squared"
gen age = sqrt(age_2)
label var age "Age"
gen pid5= 1 if k_8b==2
replace pid5= 2 if k_8b==4
replace pid5= 3 if k_8b==5
replace pid5= 4 if k_8b==3
replace pid5= 5 if k_8b==1
label var pid5 "Party ID"
label define mpid5 1"Republican"  2"Lean Republican" 3"Independant" 4"Lean Democratic" 5"Democratic"

label value pid5 mpid5


rename k_16a edu
label var edu "Education"

rename a_3b agency_size
label var agency_size "Total Agency Employees"

rename a_4b agency_budget 
gen log_agency_budget = ln(1+agency_budget) 
label var log_agency_budget "ln(Agency Budget, $2018)"

revrs e_1a 
revrs e_1b 
revrs e_1c 
revrs e_1d 

label var reve_1a "Gov."
label var reve_1b "Gov. Staff"
label var reve_1c "Legis."
label var reve_1d "Legis. Staff"



***** CODE USED TO PRODUCE APPENDIX TABLE B6 ******
ologit reve_1a i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , vce(boot, reps(1000) seed(55803)) 
est sto boot_m1

ologit reve_1b i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , vce(boot, reps(1000) seed(55803)) 
est sto boot_m2

ologit reve_1c i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , vce(boot, reps(1000) seed(55803)) 
est sto boot_m3

ologit reve_1d i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year , vce(boot, reps(1000) seed(55803)) 
est sto boot_m4

  esttab boot_m1 boot_m2 boot_m3 boot_m4 using Table_B6.rtf,starlevels( * 0.05 ** 0.01) cells((b(star fmt(3)) ci_bc[ll] ci_bc[ul] )) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Gov." "Gov. Staff" "Legis." "Legis. Staff" ) 


***** CODE USED TO PRODUCE APPENDIX TABLE B7 ******

 ologit d_15a i.intersection i.reve_1a i.reve_1b i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year, vce(boot, reps(1000) seed(55803)) 
 est sto boot_m5

 
 ologit d_16a i.intersection i.reve_1c i.reve_1d  i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year, vce(boot, reps(1000) seed(55803)) 
 est sto boot_m6

 
 ologit d_20a i.intersection i.reve_1a i.reve_1b i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year, vce(boot, reps(1000) seed(55803)) 
 est sto boot_m7

 ologit d_21a i.intersection i.reve_1c i.reve_1d  i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year, vce(boot, reps(1000) seed(55803)) 
 est sto boot_m8
 
 esttab boot_m5 boot_m6 boot_m7 boot_m8 using Table_B7.rtf,starlevels( * 0.05 ** 0.01) cells((b(star fmt(3)) ci_bc[ll] ci_bc[ul] )) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Gov. Policy" "Legis. Policy" "Gov. Regs" "Legis. Regs" ) 


***** CODE USED TO PRODUCE APPENDIX TABLE B8 ******
mlogit reve_1a i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year, base(3)  r 
est sto mlogit1

 esttab mlogit1 using Table_B8.rtf  ,starlevels(+ 0.10 * 0.05 ** 0.01) cells(b(star fmt(3)) se(par fmt(3))) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Gov." ) 


***** CODE USED TO PRODUCE APPENDIX TABLE B9 ******
mlogit reve_1b i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year, base(3)  r 
est sto mlogit2

 esttab mlogit2 using Table_B9.rtf  ,starlevels(+ 0.10 * 0.05 ** 0.01) cells(b(star fmt(3)) se(par fmt(3))) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Gov. Staff" ) 

***** CODE USED TO PRODUCE APPENDIX TABLE B10 ******
mlogit reve_1c i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year, base(3)  r 
est sto mlogit3

 esttab mlogit3 using Table_B10.rtf  ,starlevels(+ 0.10 * 0.05 ** 0.01) cells(b(star fmt(3)) se(par fmt(3))) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Legis." ) 

***** CODE USED TO PRODUCE APPENDIX TABLE B11 ******
mlogit reve_1d i.intersection i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom i.funcat13 i.state i.year, base(3)  r  
est sto mlogit4

 esttab mlogit4 using Table_B11.rtf  ,starlevels(+ 0.10 * 0.05 ** 0.01) cells(b(star fmt(3)) se(par fmt(3))) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Legis. Staff" ) 

***** CODE USED TO PRODUCE APPENDIX TABLE B12 ******
mlogit d_15a i.intersection i.reve_1a i.reve_1b i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year , base(2) r 

est sto mlogit5

 esttab mlogit5 using Table_B12.rtf  ,starlevels(+ 0.10 * 0.05 ** 0.01) cells(b(star fmt(3)) se(par fmt(3))) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Gov. Policy" ) 


***** CODE USED TO PRODUCE APPENDIX TABLE B13 ******
 mlogit d_16a i.intersection i.reve_1c  i.reve_1d i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year , base(2) r 

est sto mlogit6

 esttab mlogit6 using Table_B13.rtf  ,starlevels(+ 0.10 * 0.05 ** 0.01) cells(b(star fmt(3)) se(par fmt(3))) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Legis. Policy" ) 


***** CODE USED TO PRODUCE APPENDIX TABLE B14 ******
mlogit d_20a i.intersection i.reve_1a i.reve_1b  i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year , base(2) r  

est sto mlogit7

 esttab mlogit7 using Table_B14.rtf  ,starlevels(+ 0.10 * 0.05 ** 0.01) cells(b(star fmt(3)) se(par fmt(3))) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Gov. Regs." ) 


***** CODE USED TO PRODUCE APPENDIX TABLE B15 ******
mlogit d_21a i.intersection i.reve_1c i.reve_1d i.civil_service weekly_hours age age_2 b3.edu years_employ_state years_employ_agency years_employ_position b3.pid5 i.agency_size log_agency_budget inst6017_nom  i.funcat13 i.state i.year , base(2) r 

est sto mlogit8

 esttab mlogit8 using Table_B15.rtf  ,starlevels(+ 0.10 * 0.05 ** 0.01) cells(b(star fmt(3)) se(par fmt(3))) ///
 stats(N aic bic, labels("N. of observations" "AIC" "BIC" )) ///
 label varwidth(35) ti("") ///
 mtitles("Legis. Regs." ) 






